Recovering magnetization distributions from their noisy diffraction data 
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shot-noise and the intensity of charge scattering noise are independently varied. This work is directly relevant to 
the time-resolved imaging of magnetic dynamics using coherent and ultrafast radiation from X-ray free electron 
lasers and also to broader classes of diffractive imaging experiments which suffer noisy data, missing data or 
both. 
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I. INTRODUCTION. 

There has been a growing interest in studying and manip- 
ulating thin-film magnetic nanostructures flTJ-Q] . Besides the 
commercial applicability of such studies, experimental data 
on the formation, dynamics and stability of magnetic nanos- 
tructures will provide clues for constructing predictive models 
of magnetic materials [6], which may in turn drive the inven- 
tion of novel devices. 

A comprehensive understanding of these magnetic nanos- 
tructures involves studying extremely fast magnetic dynamics 
at high resolution. Ideally, this can be achieved by sequen- 
tially illuminating an evolving magnetic specimen using very 
short, intense pulses of coherent X-ray radiation (image of 
such in Fig. Q]). Such radiation has become available at X-ray 
free electron laser (XFEL) facilities, which can produce fem- 
tosecond pulses with upwards of 10 12 X-ray photons each. 
Despite such high intensities, pulses are typically monochro- 
matized and polarized for magnetic imaging at the expense of 
their photon flux. Furthermore, in the case of repetitive stud- 
ies on the same sample, say to study a specimen's dynamics, 
the intensity of the XFEL pulses may have to be reduced to 
prevent sample damage by energetic X-ray photons. As a re- 
sult of reducing pulse intensity, the diffraction signal from the 
weakly scattering magnetic contrast in specimens are often 
expected to be photon-shot-noise limited [4]. To make mat- 
ters worse, the magnetic signal may also be contaminated by 
strong scattering from the non-uniform charge density intrin- 
sic to magnetic specimens. 

Currently, Fourier transform holography J2l and speckle 
metrology [3] are two leading coherent X-ray techniques al- 
ready used to study magnetic nanostructures. Their effective- 
ness comes with limitations: speckle metrology is restricted 
to "fingerprinting" in reciprocal space (unable to resolve lo- 
calized dynamics of Fig. [T); Fourier transform holography 
affords direct-space imaging but it requires the crafting of a 
reference structure. 

This paper describes an alternative to Fourier transform 
holography and speckle metrology, detailing a reconstruction 



algorithm that directly images extended magnetization distri- 
butions when given only transmission diffraction data, with- 
out the need for a reference illumination. This algorithm cru- 
cially uses prior information about the ensemble of magneti- 
zation distributions to reconstruct a specific distribution. To 
demonstrate this, we had to generate an ensemble of credible 
magnetization distribution to be used as scattering sources for 
our diffractive imaging simulations (SectionlHlTi. Whereas our 
algorithm is robust when the ensemble magnetization values 
are known, it is still relevant even when given limited infor- 
mation about these values (AppendixlBl. We expect our algo- 
rithm to apply, with possibly reduced efficacy, to real magne- 
tization distributions whose ensemble properties are less well 
characterized but qualitatively similar to those in this paper. 

Our algorithm also exploits knowledge of the sample's 
direct-space support to reduce the effects of photon-shot-noise 
in the diffraction data. To make our demonstrations relevant 
to ultrafast magnetic imaging with XFEL radiation, we use 
diffraction data with the severest levels of noise in addition 
to missing data in the beamstop. We compare our diffractive 
imaging reconstruction with the performance of Fourier trans- 
form holography at comparable noise levels. 

To substantiate our methods, we include a feasibility study 
of noisy magnetic imaging when subjected to our described 
methods in the most optimistic scenario that the ensemble 
magnetization values are known. We expect this feasibility 
study to be useful in the design of future experiments. 



II. RESONANT MAGNETIC SCATTERING. 

Multilayer magnetic thin films with perpendicular magnetic 
anisotropy exhibit a notable phase comprising mag- 

netic nanostructures that can be described by a 2D coarse- 
grained magnetization distribution m(r). The magnetization 
in this phase is effectively parallel or antiparallel to the sam- 
ple's layer normal IU7I1 . In this section, we briefly discuss how 
such magnetization distributions are encoded in the diffraction 
data. 
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/ tot (r)«/ c (r) + rm(r) 
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FIG. 1: Simulated example of the type of magnetic domain motion 
that could be investigated by an XFEL source. The magnetization 
distributions above differ only within the central circle. The average 
domain is 5 pixels wide, or approximately 170 nm when related to 
the experiment data in [2] (each pixel hence measures 34 nm). We 
explain how to generate an ensemble of such similar distributions in 
the text. 



In diffractive imaging, one typically measures the sample's 
elastic photon scattering amplitude which varies across the 
sample. This scattering mechanism includes virtual transi- 
tions between core electron states and unoccupied electron 
states above the Fermi level [8]. Since these unoccupied states 
are spin-polarized due to the sample's local magnetization, the 
photon scattering amplitude depends on the sample's magne- 
tization distribution m(r). 

There are, naturally, other components of the sample's elas- 
tic scattering amplitude that are insensitive to the magnetiza- 
tion: /°(r), the Thomson contribution; / c (r), the anomalous 
charge scattering. Both of these contributions are integrated 
along the incident beam direction. Like m(r), /°(r) and 
/ c (r) are also treated as 2D distributions. 

A magnetic specimen's total elastic scattering amplitude is 
given by H as 

/ tot (r) « /°(r) + ( e/ * • e,) / c (r) - i(e/ x e*) • m(r) f M 
+ ( e/ * • m(r)) ( ei * • m(r)) f m , (1) 

where and e f are the polarization vectors of the incident 
and scattered radiation. The magnetization-sensitive scatter- 
ing amplitudes f M and f m are scaled to allow the magnetiza- 
tion to be normalized as max(|m(r)|) = 1. 

The total elastic scattering amplitude /tot( r ) of multilayer 
magnetic films can be simplified with a few experimental con- 
straints. First, since the magnetization is parallel or antiparal- 
lel to the sample's layer normal, we can replace m(r) with 
the longitudinal scalar distribution m(r). More importantly, 
the contribution to the scattering amplitude from the last term 
in eq. (Q} is suppressed if light were transmitted along this 
longitudinal direction. Second, we restrict ourselves to circu- 
larly polarized incident radiation, which is a scattering eigen- 
state of the 3rd term in eq. (fTJ. This choice, however, causes 
the diffraction patterns from magnetic and charge distributions 
to interfere [18]. Third, in the small-angle scattering limit, 
we can combine the non-magnetic scattering contributions as 
/ (r). These conditions produce a simplified total scattering 



Experimently, the magnetic scattering amplitude f M can be 
dramatically increased through resonant scattering: by tuning 
the energy of the incident photons to match those of core-level 
electron transitions in the sample (L or M edges) HI IH] . 
This enhances the scattering signature of the magnetization 
with respect to the charge distribution, which is useful since 
we are interested only in the former. 

A difference in the correlation length of the charge dis- 
tribution and that of the magnetization distribution is com- 
mon, which causes a separation in the peaks of their respec- 
tive diffracted power [9]. Potentially, one could then ignore 
the charge distribution when imaging the magnetization at a 
lower resolution. However, later sections of this paper show 
that magnetic imaging may still be difficult despite such a sep- 
aration. 

The Fraunhofer diffraction intensity from samples obeying 
eq. (O is 
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d 2 r e^ r f tot (r) 
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where cf>, the photon fluence, crucially determines the number 
of diffracted photons and hence the severity of photon-shot- 
noise. Since only the total number of diffracted photons in our 
simulations is experimentally relevant and can be varied by 
changing only </>, the absolute scale of the magnetic and charge 
scattering amplitudes in /tot(r), / and / respectively, is 
immaterial. From here on, magnetic scattering amplitudes 
and magnetization become interchangeable because they dif- 
fer only by this unimportant absolute scale. The same is true 
between charge scattering amplitudes and charge. Since it still 
serves to be consistent, we normalize the magnetic scattering 
amplitude f M = 1 in eq. (|2). The ratio f c / / , however, 
depends on the experimental specimen and the polarization 
of the incident radiation. This means, of course, that 4> is no 
longer strictly the photon fluence, but a variable to control the 
number of scattered photons. 



HI. GENERATING MAGNETIC DOMAIN PATTERNS. 

To simulate realistic magnetic imaging, we first need to 
generate magnetization distributions, or domain patterns, that 
resemble a wide and interesting variety of actual specimens. 
At a minimum, the ensemble of such domain patterns should 
conform to these experimental observations: 

1. in Fourier-space, an azimuthally symmetric diffracted 
power which peaks at a particular spatial frequency 
(compare simulated example in Fig. [2J3 to those from 
experiments in lHHSH); 

2. in direct-space, a statistical distribution on the magne- 
tization (Fig. |4]i of ferromagnetic domains with finite- 
width domain walls (Fig. |2^). 
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The clues to generating realistic domain patterns lie in the 
careful examination of the diffraction envelope shown in Fig. 

The spatial frequency dependence of this envelope reveals 
two competing effects: short-range exchange interaction that 
produces ferromagnetic domains and long-range demagnetiz- 
ing fields which in turn destabilize these domains. 

These effects are modeled by the 2D Landau-Ginzburg free 
energy density 



T(r) = A(m{v) 2 -l) +B\Vm{v)\ 2 
C 
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where A, B and C are temperature dependent positive quan- 
tities, and I = Tr/q max is a cutoff that defines the maxi- 
mum spatial frequency. Rewriting (01 in terms of the Fourier 
modes of the magnetization, m(q), we obtain in the limit 
|q| = q <C (/max the following expression: 



.F(q) = |m(q)| 2 

+ 0(\m(q)\ 



-2A + Bq 2 + 

2irC (q max /TT -q + (tt/4) g 2 /<7max) 



(5) 



Defining new constants a > and b, and rescaling g max by a 
constant, this can be rewritten in the much simplified form 

•F(q) - [(S/W ~ a) 2 + b] |m(q)| 2 + O (|m(q)| 4 ) . (6) 

The ferromagnetic instability corresponds to b — >• and the 
g-dependence of the fluctuations as this limit is approached is 
controlled by the coefficient of the term quadratic in the mag- 
netization. As a simple model for the formation of magnetic 
domains in real materials we will assume the distribution of 
fluctuations in the paramagnetic phase (b > 0), given by the 
equipartition theorem, is preserved when the system freezes 
into a particular domain pattern. The intensity in this model is 
given by 
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Our simulations will use this form for the power spectrum 
with a and b fitted to agree with experimental data in |01 • We 
use dimensionless units where the maximum spatial frequency 
9max is scaled to the value tt. 

The generation of each domain pattern begins with an ar- 
ray of random, uniformly distributed real numbers between -1 
and +1, mimicking the high-temperature magnetization distri- 
bution in the absence of external fields. On this random state 
m(r), we apply two nonphysical operations in turn: 



1 . band-pass Fourier filter using eq. (Q - 



m(q) 



m(q) 



(8) 



\/(<z/<W - 0.27) 2 - 0.015 ' 
where m(q) is the discrete Fourier transform of m(r); 




FIG. 2: (color online). Simulated magnetic domains and their 
diffraction intensities, which together illustrate the effects of charge 
scattering. Panel (a) shows domains without charge scattering and 
in (b) the logarithm of their diffraction intensities. In (c) the same 
domains with random charge distribution (A m /A c = 1) added and 
(d), the logarithm of its resultant diffraction intensities. Panel (e) 
plots the azimuthally averaged diffraction intensities from the charge 
only (thin, solid line), magnetic only (thick, solid line) and charge- 
plus-magnetization distributions (dashed line) belonging to the do- 
main pattern in (c). The inset in (e) plots the scattering amplitudes 
(ff for charge; rrn for magnetic) along the horizontal white line of 
the simulated domain pattern in (c). 



2. binary projection on m(r) 
m(r) 




if m(r) > 
if m(r) < 



(9) 
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The composition of these two operations is iterated on m(r) 
until it converges, where the values of m(r) are unchanged 
upon further iteration. Thereafter, we simulated finite domain 
wall widths by multiplying the converged distribution m(r) 
with a final low-pass Fourier filter, exp (—2.5 (q/qmax) 2 ) 11911 ■ 
This domain pattern is then normalized to max(|m(r)|) = 1. 

Different, random initial arrays result in different domain 
patterns m(r), defining an ensemble of simulated patterns. 
Whereas we generated domains with zero net magnetization, 
this recipe can be easily modified to change this net magneti- 
zation. 

This recipe for generating domain patterns is easily ex- 
tended to create perturbed versions of any domain pattern: 
we replace randomly selected circular areas in a previously 
converged source domain pattern with random numbers, then 
reapply the domain generation recipe until this perturbed pat- 
tern converges. This replacement occurs before the low-pass 
Fourier filter is applied to the source pattern. As an example, 
the pattern in Fig. QJ) is a converged perturbation of the pattern 
in Fig. QJi. 



IV. MODEL OF CHARGE SCATTERING. 

Since it is reasonable to expect the charge distribution 
/ (r) to be spatially uncorrected at the resolution of the 
resonant scattering experiments |4j,|9[], we model it as a 2D 
array of random, real numbers / c (r). Each array element 
of f c (r) represents the charge scattering amplitude averaged 
over a pixel. 

The statistics of the spatially uncorrected charge distribu- 
tions is characterized by its mean (f c (r)) and standard devi- 
ation, 

A c = ^((fC(r)-{fC(rW), (10) 

which we coin charge contrast. The angle brackets denote the 
average over each distribution. The charge contrast should be 
compared to the magnetic contrast, 

A m = V(Mr)-Mr))) 2 ) ■ (11) 

The diffraction intensity in eq. (01 does not distinguish be- 
tween charge and magnetic scattering, so any reconstruction 
can only determine their sum (see eq. (ffj). Since we are 
interested only in recovering the magnetization distribution, 
the intrinsic charge distribution will contribute an inextrica- 
ble scattering noise, characterized only by the signal-to-noise 
ratio A m /A c . When A m ks A c , it becomes visually impos- 
sible to differentiate between these distributions even if their 
sum were correctly reconstructed (compare Fig. [2^ and|2};). 

In contrast, the mean charge scattering amplitude (f c (r)), 
as later sections will show, is an immaterial constant to the 
reconstruction of m(r). Nevertheless, to be consistent, we 
fixed (f c (r)) = 1.33 using experimental data from yfl. 
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FIG. 3: (color online). Noisy photon data from simulated diffraction 
experiments. Panel (a) shows the diffraction data from the domain 
pattern in Fig. [9^, with signal-to-noise of point B in Fig. [8] The 
intensities in the dashed-line box (lower right) of panel (a) are inver- 
sion symmetric to those in the solid-line box (upper left); larger inset 
is a magnified view of photon data in the solid-line box; the central 
black disk is the beamstop. Panel (b) shows the azimuthally averaged 
photon counts in panel (a). 



V. DIFFRACTIVE IMAGING AS CONSTRAINT 
SATISFACTION. 

We can interpret the diffractive imaging experiments in the 
language of constraint-satisfaction problems. Essentially, the 
goal of diffractive imaging is to recover the true magnetiza- 
tion distribution subject to two constraints: its measured noisy 
diffraction data (Fourier constraint) and the assumed statistics 
on its expected magnetization (direct-space constraint). This 
section discusses how we generated and characterized these 
two constraints. 



A. Fourier constraint. 

The Fourier constraint requires that the diffraction intensi- 
ties of the true magnetization distribution, which we wish to 
recover, be statistically compatible with the measured photon 
data, mindful that the data includes intrinsic charge scattering. 

To simulate the diffraction data, we first added each pair of 
randomly-generated magnetization and charge distributions, 



5 



m(n) 0. - 
-0.5 
-1. _ 




0% 50% 100% 

position along sorted list 

FIG. 4: Sorted-value magnetization constraint. Sorted list of nor- 
malized magnetization fh(n) when averaged over many simulated 
domain patterns (curve) and those from one pattern (dots). 



m(r) and / (r) respectively. We confined this total scatter- 
ing amplitude to a circular support S (Fig. [2^, for example). 
Its continuous intensity distribution was scaled by <f> to give 
the desired average number of scattered photons, then Pois- 
son sampled to simulate photon-shot-noise. Following this, 
we averaged each data with its Friedel-symmetry counterpart 
to make it consistent with the real-valued direct-space con- 
trast. Finally, a beamstop was applied to this symmetrized 
data, thus removing photon counts that would be contami- 
nated by intense, unscattered radiation in actual experiments. 
The size of the beamstop was adjusted such that the remaining 
photon counts span two orders of magnitude (example photon 
data in Fig. [3]). Naturally, the unmeasured Fourier amplitudes 
at spatial frequencies within the beamstop are unconstrained 
in our reconstruction algorithm. 



B. Direct-space constraint. 

When discussing the direct-space constraints on the mag- 
netization it is convenient to introduce the sorted-value rep- 
resentation m(n), where m(l) is the smallest magnetization 
among the pixels within the support, m(2) is the next smallest, 
etc. and m(Ns) is the largest magnetization value. Within the 
ensemble of random domain patterns produced by the same 
magnetic material (and identical external parameters) the plots 
of the functions m(n), with n ranging from 1 to Ng, should be 
nearly the same. Figure [4] compares m(n) for one simulated 
domain pattern with the averaged m(n) over many patterns. 
The structure of m(n) is mainly a function of two lengths: the 
width of domains and the width of domain walls; materials 
with very thin domain walls will have a more step-like m(n) 
& 

The simplicity of the function m(n) serves as a powerful 
constraint for the domain reconstructions. We will use the 
notation rh(n) for the sorted-value magnetization constraint. 
In most of our simulations we will assume m(n) is known. At 
zero net magnetization we have the symmetry 



m(n) = -m(N s + 1 - n) . (12) 

Appendix iBl briefly describes reconstructions that only use 
this property and boundedness of m(n), rather than an explicit 
constraint function rh(n). 

The knowledge of m(n), which includes information about 
the size and shape of the support S, constitute the direct-space 
constraint in our reconstructions. 



C. Noisiness of constraints. 

To prepare for systematic studies of reconstruction feasibil- 
ity, we classify our diffractive imaging simulations using con- 
venient signal-to-noise parameters. One such consideration is 
the photon-shot-noise in diffraction data. 

Photon-shot-noise is related to the average number of scat- 
tered photons per pixel ^ tot , regardless of whether it came 
from the charge or magnetization distribution (refer to eq. 
(f3]l). Increasing /i tot ought to improve the chances of recon- 
structing the total scattering distributions. However, using /i to t 
as a signal-to-noise parameter is too optimistic since we are 
only interested in recovering the magnetization distribution 
112111 . Consequently, one must still isolate the magnetization 
distribution from the total scattering distributions, even if the 
latter is correctly determined (i.e. to extract the magnetization 
Fig. [2^ from only Fig. [2}:). 

To appropriately characterize the noisiness of the photon 
data to our goal, instead of the total scattered power /i to t, we 
use n m : the average number of photons scattered due to the 
magnetization in each pixel within the support S. In exper- 
iments, fi m can be estimated directly from magnetic elastic 
scattering amplitude f M , photon flux and exposure time of 
diffraction measurement. In our simulations, \i m is computed 
as 

N S ^ 

where m(q) is the discrete Fourier transform of the magneti- 
zation distribution m(r); N$ is the number of support pixels; 
(\> is the same scalar in eq. (01 which we vary to give the de- 
sired total number of scattered photons; f M is again set to 
unity inconsequentially. The product ^L m Ns corresponds to 
the total number of photons scattered per pulse in the absence 
of charge scattering. 

The other noise consideration comes from charge scatter- 
ing. We assume that the specimen's random charge distri- 
bution is unknown, which results in a harder reconstruction 
problem. As a consequence, the model magnetizations in Fig. 
|4]will not agree with those in the total scattering amplitudes 
of eq. ©, which includes the charge distribution. Essentially, 
this makes our direct-space constraint noisy l22ll . Experimen- 
tal measurement of the charge distribution would certainly re- 
duce this noise and simplify the reconstruction. 
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VI. RECONSTRUCTION ALGORITHM. 
A. Modifying the difference map. 

Seeking the true magnetization distribution is equivalent to 
finding the intersection of the Fourier and direct-space con- 
straint sets. Such intersections, or solutions, can be discov- 
ered using an iterative constraint-satisfaction algorithm: the 
difference map [10], which uses simple projections to these 
two constraints (Pd, projection to direct-space constraint; Pp, 
projection to Fourier constraint). 

The difference map algorithm accelerates the discovery of 
a solution, primarily by reducing the dimension of the search 
space [10]. It is also particularly efficient in extricating the 
iterate from near intersections (false solutions) to prevent the 
search from stalling. However, the difference map algorithm 
was optimized for noiseless constraints sets with true intersec- 
tions JH. 

Unfortunately, photon and charge scattering noise distorts 
our measurement of the true Fourier constraint, demoting its 
intersections with the direct-space constraint to near intersec- 
tions, from which iterates are jettisoned. This prohibits the 
search from reaching the true magnetization distributions en- 
coded in these near intersections. 

To increase its reconstruction success rate, the difference 
map was modified to improve the stability of the iterate around 
a near intersection. This is accomplished by an intermediate 
step to the iteration rn n — > m' n — > to„ +1 (where the iteration 
number n = 0,1, 2,. . .), which keeps the iterate close to the 
Fourier constraint 02311 : 

m' n = am n + (1 - a)P F (m n ) , 
e n = Pp(2P D (m' n )-m' n )-P D (m' n ) , 
m n+ i = m' n + e n , (14) 

with Pjj and Pp as the direct-space and Fourier constraint 
projections respectively and a as the map's modification pa- 
rameter. The update on the iterate is denoted e n , so that it may 
be referenced concisely in later paragraphs. 

In our reconstructions we chose a — 0, which substan- 
tially improves the iterate's stability (Fig. [5]l while reducing 
the number of computations in the first step of the algorithm. 
With a = 1, eq. ( TT~4T > reduces to an instance of the original 
difference map. AppendixlAldiscusses how eq. (fl4l is similar 
to the RAAR algorithm in ifiai. 

The modified difference map is iteratively applied to a ran- 
dom, initial magnetization distribution mo- The norm of the 
map's update | \e n \ \, which we term the error metric, measures 
the average change of the iterate during the search. When the 
error metric drastically declines, it indicates that the differ- 
ence map updates have experienced a dynamic transition and 
the search has likely converged (Fig. |6). Because of the inher- 
ent noise in the constraints, the error metric will never vanish 
as it would, had an intersection of the two constraints been 
found in the noiseless case. When a noticeable transition in 
1 1 e ra 1 1 occurs and is stable, we harvest the candidate solution of 
the magnetization distribution, Pd (Pp ( 2Pp>(m' n ) — m' n )). 
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iteration number, n 

FIG. 5: Stability of modified difference map, eq. J 1 4b . around a so- 
lution. We estimate this stability by tracking how quickly the iterate 
leaves the solution because of the noisy constraints. Starting from the 
solution magnetization distribution, the normalized magnitude of the 
iterate's updates (||£n||/-/V^ 2 ) remains low for a < 0.5, indicating 
iterate stability around a noisy solution. 
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FIG. 6: Signature of a successful reconstruction. The normalized 
error metric (||£n||/A r g / ' 2 ) suffers a noticeable transition at iteration 
number n « 50000 during the successful reconstruction of the dis- 
tribution in Fig. [9j). 



The correctness of this candidate solution is tested when com- 
pared against other candidate solutions from different, ran- 
dom, initial iterates too. Consistent recovery of nearly identi- 
cal candidate solutions, up to an overall multiplicative sign or 
spatial inversion, from random restarts asserts their credibil- 
ity as the true magnetization distribution. One can smooth out 
the fluctuations between the candidate solutions by averaging 
them. 

In searches using the noisiest photon data, the error metric 
\\e n 1 1 will never show a clear transition. In such cases, recov- 
ering the true magnetization distribution is plainly impossible. 
Nonetheless, we can still evaluate the search results, however 
wrong they may be. From eq. (fT4b . notice that \ \e n \ \ also mea- 
sures the distance between two points on the two constraints: 
Pp ( 2Pp,(m' n ) — m' n ) and Pp)(m' n ). Hence the minimum 
\\e n \ \ during a search signals the nearest distance between the 
two constraints — the best alternative to discovering an in- 
tersection. Unlike more robust candidate solutions with less 
noisy data, these faux solutions are never repeated with ran- 
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photon data healed photon data true intensities 



FIG. 7: (color online). Speckle-healing by applying an autocorrela- 
tion support constraint to the photon data during reconstruction. We 
show this healing of a magnified section of the photon data in Fig. [3] 



dom restarts. 



indirectly constrained by the diffraction data and after nu- 
merous iterations the difference map iterate proposes prelimi- 
nary intensities for them. We replace the missing photon data 
with these preliminary intensities before applying the speckle- 
healing operation. In our reconstructions, the photon data was 
healed this way every 1000 iterations, which were then used 
to constrain iterations until the next healing. When the recon- 
struction converges under this adiabatic healing process, we 
are assured that it is still compatible with the photon data. 

With this adiabatic speckle-healing procedure in effect, the 
projection to the Fourier constraint, Pp(m), comprises the 
following operations on the iterate's Fourier transform m(q): 

1. set the amplitudes of m(q) outside the beamstop to the 
square root of the speckle-healed photon data, while re- 
taining the phases of m(q); 



B. Projection to direct-space constraint. 2. the m(q) values within the beamstop are unchanged. 



The projection to the direct-space constraint, Pn(m), com- 
prises the following operations on m: 

1. set all values of m(r), for r outside S, equal to zero; 

2. m(n) replaced by the magnetizations m(n) shifted and 
scaled to have the same mean and variance as m(n) be- 
fore the projection. 

Step |2] allows the mean scattering amplitude (/ c (r)) and 
the magnetic contrast A m to be indirectly constrained by the 
diffraction data 12411 . 

In actual experiments where the magnetization constraint 
function fh(n) is not readily available or simulated, one could 
instead project m(n) to a class of parametrized magnetization 
functions, where the projection determines the best parameter. 
When even this is impossible, imposing only the key features 
of fh(n) on m(n) may be sufficient (AppendixlBl. 



C. Projection to Fourier constraint. 

Before discussing the projection to the Fourier constraint, 
Pf(w), we describe a modification to the diffraction data 
which lowers the photon-shot-noise using the direct-space 
constraint. If the scattering distribution is contained within 
a direct-space support S, the Fourier transform of the diffrac- 
tion intensities — or the autocorrelation of the direct-space 
scattering distribution — should be contained within the auto- 
correlation support Sa- 

We can lower the noise in the diffraction data using our 
knowledge of the support, hence constraining the photon data 
to have the expected speckles. We did so by applying an au- 
tocorrelation support constraint to the Fourier transform of 
the photon data — setting all values outside Sa in the data's 
Fourier transform to zero. Empirically, this speckle-healing 
operation increases the R-factor between the processed pho- 
ton data and the true intensities (see Fig. [7J. 

However, the missing data within the beamstop may con- 
fuse speckle-healing. These central Fourier amplitudes are 



VII. FEASIBILITY. 



A. Difference map reconstructions. 

Unlike an actual experiment, the true magnetization distri- 
butions are known in our simulated experiments. This allows 
us to directly compare the reconstructions m rec with the true 
distribution m true within the support S via the following devi- 
ation measure: 

s = jEK ec ( r )- wtme ( r )) 2 < 15 > 



res 



u r 



m(r) 



(Eres™( r )' 



,1/2 



(16) 



The deviation i5 is proportional to the square of the distance 
between u lec and it 1 ™ 6 , which are the respective distributions 
normalized as unit vectors. Allowing for an overall sign in the 
reconstructed magnetization, deviations lie within the range 
< S < 1. 

In our simulations, reconstructions with 5 < 0.2 were con- 
sistently recovered from random restarts. In actual experi- 
ments, only the consistency test is available to evaluate the 
reliability of the reconstructions. We deem such repeatable 
reconstructions to be successful. 

We systematically studied the performance of our recon- 
struction algorithm when we varied the two signal-to-noise 
parameters: relative magnetic contrast A m /A c and the aver- 
age number of photons scattered from the magnetic distribu- 
tion fi m . As Fig. [Syndicates, increasing A m /A c and p, m im- 
proves the accuracy of the reconstructions. However the same 
figure shows that the effects of suppressing either A m /A c or 
p, m are qualitatively different — lowering one variety of noise 
can not compensate for the reconstruction errors caused by the 
other. 

Magnetization distributions shown in Fig. [9] are routinely 
reconstructed with low deviations 6 < 0.1 given remark- 
ably noisy data typical of Fig. [3j relative magnetic contrast 
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FIG. 8: (color online). Feasibility of diffractive imaging at various 
noise levels. The top panel shows a contour plot of the reconstruc- 
tion deviation 8 (in eq. d!5t ) generated from many simulated re- 
constructions of the pure magnetization distribution similar to (a) 
as the signal-to-noise parameters (A m /A c and fj, m ) were indepen- 
dently varied. Panels A, B and C show reconstructions of (a) subject 
to corresponding signal-to-noise parameters marked in the top panel. 



A m /A c = 5 and average scattered photons due to the mag- 
netization in each pixel /i m = 3. The deviation of reconstruc- 
tions from the true domain pattern at various noise levels is 
numerically computed in Fig. [8] and appears to be indepen- 
dent of the support size at a constant domain resolution (the 
reconstructions in Fig. [8)3 and Fig. [9J) suffered comparable 
noise levels). 

Reconstructions with the unmodified difference map, a 
1 in eq. ( fT4l . do not converge within the range of noise pa- 
rameters in Fig. [8] neither in the sense of achieving a dy- 




FIG. 9: In panel (b), a modified difference map reconstruction of a 
simulated magnetization distribution, panel (a), with signal-to-noise 
parameters corresponding to point B in Fig. [8] The diffraction data 
used for this reconstruction is shown in Fig. [3] If the magnetic do- 
mains are 170 nm wide, then having only three photons scattered 
by the magnetic contrast within each 34 nm pixel was sufficient to 
reconstruct the pattern in the bottom panel. 



namic transition in the error metric (Fig. |6]l nor repeatability 
given random restarts. We witness this lesser performance 
even with reconstructions using the modified difference map 
when we omit either the sorted-value magnetization constraint 
in the direct-space projection or the speckle-healing procedure 
or both. 
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FIG. 10: Simulated reconstructions with Fourier transform hologra- 
phy (FTH). The panels (b), (c), (d) show how FTH require more scat- 
tered photons for acceptably low reconstruction deviations. Panel (e) 
is a FTH reconstruction of Fig. [9^. Panels (d) and (e) together show 
that reconstructions worsen with increasing support size. The rela- 
tive magnetic contrast was A m /A c = 10 in these reconstructions. 



B. Comparison with Fourier transform holography. 

To provide perspective, we compared our reconstructions 
with those from simulated Fourier transform holography 
(FTH) in Fig. |T0j In FTH, the domain pattern is obtained 
directly from its cross-correlation with an aptly machined ref- 
erence pinhole |2y. This cross-correlation is obtained from a 
simple Fourier transform of the measured diffraction intensi- 
ties without the need for phase retrieval. 

To make the comparison more compelling, we provided our 



FTH simulations with the following advantages over the non- 
holographic method: 

1. noisy diffraction signal within the beamstop region was 
provided; 

2. single-pixel reference pinhole for highest possible re- 
construction resolution (pinhole diameter roughly 34 
nm if magnetic domains are 170 nm wide) whereas the 
pinhole in [2] which had an effective X-ray transmis- 
sion diameter of approximately 100 nm. 

In our simulated FTH reconstructions, the ratio of the number 
of photons scattered by the magnetic contrast in each support 
pixel to the number which pass through each pixel of the ref- 
erence pinhole is 1:50, as estimated from |2|]. 

At the low signal-to-noise levels of Fig. [8j the low de- 
viation reconstructions using our proposed non-holographic 
diffractive imaging technique are out of the reach of our im- 
plementation of FTH (compare Figs. [51 to ITOb. Fig. [TOlalso 
illustrates that our FTH reconstructions worsen with increas- 
ing support size because the photon fluence through the pin- 
hole does not increase commensurate ly [25]. This reflects the 
typical situation in microscopy that higher resolution neces- 
sitates a smaller field of view. Non-holographic diffractive 
imaging does not suffer this size dependency since only the 
noise per support pixel is important. While non-holographic 
diffractive imaging does not need the experimental fabrication 
of a small reference object and can use the beam's spatial co- 
herence more efficiently via a tighter X-ray focus, it requires 
accurate knowledge of the support [13]. These differences 
between the techniques make the non-holographic phase re- 
trieval approach demonstrated here of particular interest for 
situations where the signal is too noisy for successful FTH of 
extended magnetization distributions. 

VHL CONCLUSIONS. 

Ultrafast imaging of magnetic nanostructures is presumably 
possible within the noise limits predicted by Fig. [8] This, 
of course, is valid only in the absence of other varieties of 
noise. Our study is limited to magnetic imaging without prior 
measurement of the random charge distribution. We speculate 
that the reconstruction noise limits would improve if the spec- 
imen's charge distribution, which may fluctuate, were avail- 
able. 

Certainly, imposing ensemble properties of the domain pat- 
terns in our reconstruction algorithm allows magnetic imag- 
ing with remarkably noisy data. Although our reconstruc- 
tions use the sorted-value magnetization constraint, an ap- 
proximate knowledge of this constraint may be satisfactory 
(see Appendix iBt. Despite restricting our simulations to a 
small ensemble of domain patterns, the methods we used to 
reconstruct these patterns should be valid for imaging a larger 
ensemble of ferromagnetic contrast that differ only qualita- 
tively from ours. 

Our success with the modified difference map, eq. ( TBI ), 
suggest its relevance to constraint-satisfaction problems that 
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suffer from imprecise or noisy constraints. Similarly, the 
speckle-healing procedure in this paper is pertinent to recov- 
ering missing global information common in diffractive imag- 
ing. 
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Appendix A: Similarity of modified difference map to relaxed 
averaged alternating reflections algorithm. 

The modified difference map in eq. ( fT4b resembles the re- 
laxed averaged alternating reflections algorithm (RAAR) used 
in iterative phase retrieval [12]. Like the modified difference 
map, RAAR was designed to stabilize iterates in the domain 
of attraction of a solution given noisy diffraction data. To see 
their resemblance, we combine the last 2 lines of eq. (fl4] > as a 
single operation D: 

m' n = am n + (1 - a)P F (m n ) , (Al) 
m„+i = D(m' n ) . (A2) 

The first step in the next iteration would be 

m' n+1 = am n+ i + (1 - a)P F (m n+ i) (A3) 
= aD(m' n ) + (1 - a)P F (D(m' n )) , (A4) 

which is similar in structure to the RAAR update: 

m n+ i = a D(m n ) + (1 - a)P F (m n ) . (A5) 
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FIG. 11: Normalized error metric for a successful reconstruction us- 
ing the weaker symmetry-and-boundedness value constraint. The dy- 
namic transition in the error metric is still visible, though smoothed 
out because of this less restrictive value constraint. 
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Appendix B: Symmetry-and-boundedness constraint. 

For cases when the ensemble's sorted-value magnetization 
m(n) is unavailable, it could be replaced with a less restrictive 
value constraint. One such replacement is the magnetization's 
expected symmetry in eq. (I12t . This occurs in magnetic sam- 
ples of zero net magnetization in the absence of external fields. 
In addition to this symmetry constraint, the magnetization val- 
ues must bounded by — 1 < m(n) < 1 given our normal- 
ization. Magnetization symmetry-and-boundedness together 
constitute a weaker direct-space value constraint; it is weaker 
because it includes magnetization functions besides the true 
one. 

To test its effectiveness, this weaker constraint was used to 
reconstruct a domain pattern instead of using its true model 
magnetization m(n). This mimics the experimental scenario 



FIG. 12: A reconstruction using the weaker symmetry-and- 
boundedness direct-space constraint. In panel (a), the ordered re- 
constructed magnetizations (dark curve) approximates the true val- 
ues (light curve). Panel (c) shows such a reconstruction of the source 
domain pattern in panel (b), with fi m = 4 and A m /A c = 6 (top 
right corner of contour plot in Fig. [8j. Compare panel (c) of this 
Figure with Fig. [8ji, which uses the sorted-value magnetization 
constraint and reconstructs the pattern with lower deviation while 
needing slightly fewer photons and tolerating more charge scattering 
noise. 



when one attempts reconstruction even when m(n) is not 
completely known beyond its symmetry and boundedness. 
Projecting to this constraint modifies the direct-space projec- 
tion of Section lVIBl 
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m(n) -> < _ 



sgn(m(n)) 
m(n) 



if |m(n) > 1 
otherwise, 



where 



m(n) = — (m(n) — m(N,s + 1 — n)) . 



(Bl) 



(B2) 



Practical reconstructions using this weaker direct-space 
constraint (Fig. fT2T> require the data to have slightly better 
signal-to-noise than reconstructions using the ensemble's list 
of model magnetizations (Fig. [8}. This is because the weaker 
constraint permits model magnetizations different from the 
true one and is hence a lesser guide during our search for the 
solution. 



A weaker constraint also causes smoothing of the transi- 
tion in the error metric during a successful reconstruction due 
to the relaxation of the reconstructed magnetizations towards 
the true magnetization constraint function (Fig. ITTb . With 
more photons and less charge scattering noise, the magne- 
tizations reconstructed using the symmetry-and-boundedness 
constraint (Fig. [T2a ) become closer to the true magnetization 
function. This suggests that one could obtain accurate magne- 
tization functions from low-noise diffractive imaging experi- 
ments to be used as constraints for noisier ultrafast imaging. 



We note that the symmetry-and-boundedness constraint has 
a crucial difference from the sorted-value magnetization con- 
straint: the former does not explicitly reject charge scattering 
as a source of noise while the latter does. 
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